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Abstract 

Using about 23M BB events collected in 1999-2000 with the BABAR detector, we report on the 
decays B — > D^*^D^*^ K and B^ D*^D*^ . The branching fractions of the low background decay 
modes B° D*-D'^*^'^K+ are determined to be B{B^ D*-D^K+) = (2.8 ± 0.7 ± 0.5) x lO'^ 
and B{B^ ^ D*-D*^K+) = (6.8 ± 1.7 ± 1.7) x lO'^, where the first error quoted is statistical and 
the second systematic. Observation of a significant number of candidates in the color-suppressed 
decay mode -B+ D*+D*^K+ is reported with a branching fraction B{B^ D*^D*~K+) = 
(3.4 ± 1.6 ± 0.9) X 10^3. Decays of the type B D^*W^*^ can be used to provide a measurement 
of the parameter sin2/3 of the Unitarity Triangle. For this decay mode we measure a branching 
fraction of B{B^ —>■ D*^D*^) = (8.0 it 1.6 it 1.2) x 10^^. All results presented here are preliminary. 
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1 Introduction 



Decays of B mesons that include a charmed and an anti-charmed meson are expected to occur 
through the 6 to c quark transitions b — > cW^ , where the materiahzes as a cs pair. These 
transitions are responsible for most of the Dg production in B decays. Dg production has been 
thoroughly studied in experiments running at the T(4S') resonance [|l|, ^, ^. The inclusive rate 
for Dg production in B decays was recently measured by BABAR, where a preliminary branching 
fraction I3{B DgX) = (10.93 ± 0.19 ± 0.58 ± 2.73) x 10"^ the third error being due to the 

Df <I>7r^ branching fraction uncertainty, is reported [Q. 

— +(*) 
Until 1994, it was believed that the cs quarks would hadronize dominantly as Dg mesons. 

Therefore, the branching fraction h — > ccs was computed from the inclusive B ^ Dg X, B ^ 

(cc) X and B ^ ScX branching fractions, leading to B{b — > ccs) = 15.8 ±2.8% [^]. Theoretical 

calculations are unable to simultaneously describe this low branching fraction and the semileptonic 

branching fraction of the B meson Q. It has been conjectured that B{b ccs) is in fact 

larger and that decays B — > DDK (X) (where D can be either a D^ or a D^) could contribute 

significantly. This might also include possible decays to orbitally-excited Dg mesons, B D' 'd:\ 
followed by D** ^D'^*^ K. 

In BABAR, the high statistics available allow comprehensive investigations to be made of the 
b ccs transitions. In the analysis described in this paper, we present evidence for the decays B — > 
D(*W^*^K°s and B D^*W'-*^K, using events in which both D's are completely reconstructed. 
After describing the data sample and the event selection, we show the D^*^D^*^ K signals for the 
sum of all B submodes. The branching fractions for some of the cleanest modes, such as B^ — > 
D*~ D^*^^ Q, are computed. Observation of several candidates in the color-suppressed decay 
mode B^ D*^ D*^ is also reported. 

One of the most important goals of the BABAR experiment is to precisely measure the angles of 
the Unitarity Triangle. While the decay B^ — > J/ip can be used to measure sin2/3, the Standard 
Model predicts that the time-dependent CP violating asymmetries in the decays B" jji*)- 
can also be used to measure the same quantity. An independent measurement of sin2/? in these 
modes would therefore provide a consistency test of CP-violation in the Standard Model. The 
vector-vector decay B" — > D*~^D*~ is not, however, a pure CP eigenstate. A sizeable dilution of 
the measured asymmetry may be produced by a non-negligible B-wave CP-odd component. The 
dilution can, in principle, be completely removed by a time-dependent angular analysis of the decay 
products p. 

A detailed description of the BABAR detector can be found elsewhere |^. 

2 Analysis 

Since the B mesons are produced via e'^e~ — s- T(4S) BB, the energy of the B in the T{4S) 
frame is given by the beam energy E^^^^, which is measured much more precisely than the energy 
of the B candidate. Therefore, to isolate the B meson signal, we use two kinematic variables: AE, 
the difference between the reconstructed energy of the B candidate and the beam energy in the 

center of mass frame, and m-Es, the beam energy substituted mass, defined as tues = \/^^am~~P^ 
where p*^ is the momentum of the reconstructed B in the T(4S') frame. Signal events will have AE 
close to and mss close to the B meson mass, 5.729 GeV/c^. 

^Charge-conjugate reactions are implied throughout this note. 



2.1 B D^*W^*'K decays 

The B° and B+ mesons are reconstructed in a sample of multihadron events for all possible DDK 
modes, namely 5° ^ D^*^- D^*^^ K+ , D^*)- D^*)+ , ;d^*)o^{*)o^o ^+ ^ D^*^° D^*^+ , 
d'^*^^D^*)OK+, D^*^^ D^*^~ K+ . The D° and mesons are reconstructed in the modes L>° 
K^TT^, K^TT^TT^, K^TT^TT^TT^ and D^ — > i^~7r+7r^". Charged kaon identification, with information 
from the Cherenkov angle in the DIRC and from dE/dx measurements in the drift chamber and 
in the vertex detector, is required for most D decay modes, as well as for the from B's. 

D* candidates are reconstructed in the modes D*~^ _D°7r+, D*^ D^ir^ and D*^ -0^7, 
by combining a D^ candidate with a vr", n^, or photon. Partial reconstruction of D*^'s (no vr*^ or 
7 reconstruction) is also used in the D*^ D*^K^ mode, as explained below. B candidates 

are reconstructed from the D^*\ D^*^ and K candidates. A mass constraint is applied to all the 
intermediate particles {D^, D~ , K^). When several candidates are selected per event in a specific 
B submode {e.g. B^ D^D^K'^), a value, taking into account the difference between the 
measured and the PDG values of the D masses and of the AM (for Z)*'s) is constructed and only 
the candidate with the lowest value is kept for the given submode. 



2.1.1 Evidence for signal in the sum of all B submodes 

We present here the distributions obtained by summing all possible B D^*^ D^*^ K decay channels, 
for neutral and charged B decays respectively. Since multiple candidates are removed only submode 
by submode the same event can appear several times in distributions obtained by summing over 
all modes. An event will appear in the peak near OMeV when reconstructed correctly, in the peak 
around -160 MeV if it is a D*DK {D*D*K) decay reconstructed as DDK {D*DK), and near the 
peak around +160 MeV if it is a DDK {D*DK) decay reconstructed as D*DK {D*D*K). About 
120 B^'s and 180 B^ decays have been reconstructed. The mss distributions (Figs. 1 and 2) 
contain only events with \AE\ < 24 MeV. The mss spectrum of B'^ and B^ events can be fitted 
by the sum of a background shape and a Gaussian function used to extract the number of signal 



events. The background is empirically described by the ARGUS function |10|. 



2.1.2 Measurement of exclusive branching fractions 

In this section, we present measurements of branching fractions for the three decay channels B^ —^ 
D*-D^K+, B° D*-D*^K^ and B+ D*-D*+K+. Several candidates are also observed in 
the CP eigenstate B'^ D^*')^ D^*^^ Kg but without extracting branching fractions. 

The selection efficiencies for each mode were obtained from detailed Monte Carlo simulation 
11 1 . Typical efficiencies range from 10%, for B^ D*-D^K+ with both L»°'s decaying to Kn, to 



less than 1%, for B^ D*^ D*^ K^ with D^^s decaying to Ktttt^ or XSvr. Systematic errors ac- 
count for the uncertainties on tracking and tt^ reconstruction efficiencies, K identification efficiency, 
D and B vertexing requirements, efficiency of the requirement on IS.E used to define the signal box, 
efficiency of the D mass requirement, uncertainty on the background shape, uncertainties on the 
D and D* branching fractions, uncertainties on the selection efficiencies arising from Monte Carlo 
statistics, and uncertainty on the number of produced BB events in the data sample. For the decay 
B~^ — > D*~D*~^K^a fit is performed with the sum of a Gaussian function for the signal and an 
ARGUS function for the background. The number of signal events is 8.2 it 3.5 and the number of 
background events given by the ARGUS function is 1.7. The probability that the signal arises from 
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Figure 1: tties distribution for the sum of all neutral modes 



a background fluctuation is 1.4 x 10 ^ (> 5a). The corresponding preliminary branching fraction 
is measured to be B{B+ D*- D*+ K+) = (3.4 ± 1.6 ± 0.9) x 10"^ 

In the analysis of the decay S° D*-D^*^'^K+ we require that either the or the D decays 
to Kit and we do not explicitly reconstruct the vr'^ or the photon from D*^ D^tt^ or D^^. Events 
containing D*~ D^K~^ decays are selected by requiring | Ai^j < 25 MeV. The number of signal 

events is found to be 29.6ib7.2. After correcting for the selection efficiencies and for the intermediate 

and D*^ branching fractions [12|, the preliminary branching fraction for D*^ is 

found to be B{B^ D*-D^K+) = (2.8 ± 0.7 ± 0.5) x lO'^. Events containing B^ D*-D*^K+ 
decays are selected by requiring |A£^ + 154| < 60 MeV. The average position and width of /S.E for 
B^ — s- D*~ D*^ is found to be in good agreement with expectations from B^ D*^ D*^ 
signal Monte Carlo studies. The number of signal events found is 80.2 it 15.3. To extract the 
B^ D*^ D*^ branching fraction, the contamination from decays B^ D*^ D*^ , where 
the TT~^ from the D*~^ is not reconstructed, needs to be subtracted. This contribution has been 
estimated by performing the B^ — > preanalysis on B~^ D*~ D*'^ K'^ signal Monte Carlo, 

assuming the B~^ — > D*^ D*^ K'^ branching fraction determined above. The B^ D*^D*^K^ 
background contribution is estimated to be 20.6 it 9.7 events. After subtracting this contribution, 
the preliminary B^ D*~D*^K~^ branching fraction is determined to be: i3(i?° D*^ D*^K~^) = 
(6.8 ± 1.7 ± 1.7) X 10-3. 

As pointed out in |13|, the channel B^ — > D*^ D*~ K'^ is a CP eigenstate that could be used 
for sin 2/? measurements. For decay channels B^ — > D^*^^ D^*^^ Kg involving D^^s, at least one 
decay Kit was required. The fitted number of signal events is 10.1 =b 3.7 with an estimated 




Figure 2: uies distribution for the sum of charged neutral modes 



background of 3.4 events. The probability that the signal is a fluctuation of the background is 
1.4 X 10"^ (> 5a). Most of the signal is due to the channels D*+D~K^ (4.7 ± 2.2 events 

with a background of 1 event) and 5° D*^D*'Kg (4.8 ± 2.2 events with a background of 0.3 
event) . 



2.2 S° D*+D*- decay 

mesons are exclusively reconstructed by combining two charged D* candidates reconstructed 
in a number of D* and D decay modes. The decay modes of the and D"*" used in this analysis 
were selected by an optimization of S^/{S + B) based on Monte Carlo simulations, where S is the 
expected number of signal events and B is the expected number of background events. and 
I?"*" candidates are reconstructed in the modes K~Tr~^, K~Tr~^Tr^, '7r"'"7r"'"7r , 

A-Ott+tt-, D+ K-1T+TT+, D+ KliT+ and D+ K-K+tt+. The D*+ mesons are 
reconstructed in their decays D*^ — > D^tt"*" and D*^ D+vr^. We include for this analysis the 
decay combinations D*+ D*~ decaying to (L)° 7r+, vr") or {D^ tt"*", D~ 7r°), but not {D^ 
TT^, D~ TT^) due to the smaller branching fraction and larger expected backgrounds. and 
candidates are subjected to a mass-constraint fit and then combined with soft pion candidates. 
Charged kaon candidates are required to be inconsistent with the pion hypothesis, as inferred from 
the Chcrcnkov ring measured by the DIRC and the dE/dx as measured by the SVT and DCH. To 
select B^ candidates with well reconstructed D* and D mesons, we construct a that includes all 
PDG and measured D* and D masses. In events with more than one B° candidate, we chose the 



candidate with the lowest value of x^- 

The signal region in the AE vs. m-ES plane is defined to be jA£^| < 25MeV and 5.273 < 
mES < 5.285 GeV/c^. The width of this region corresponds to approximately ±2.5(7 in both AE 
and TTT-ES- In addition only events with a < 20 are selected. To estimate the contribution from 
background in the signal region, we define a sideband in the AE vs. mss plane as | A£^| < 200 MeV, 
5.20 < mES < 5.26GeV/c2 and 50 < \AE\ < 200 MeV, 5.26 < mss < 5.29GeV/c2. We parameterize 
the shape of the background in the AE vs. tti-es plane as the product of an ARGUS function in 
mES and a first-order polynomial in AE . Based on this parameterization we estimate that the 
ratio of the number of background events in the signal region to the number in the sideband region 
is (1.72ib0.10) X 10~^. Figure 3 shows the events in the AE vs. tues plane after all selection criteria 
have been applied. The small box in the figure indicates the signal region defined above, and the 
sideband is the entire plane excluding the region bounded by the larger box outside the signal region. 
There are a total of 38 events located in the signal region, with 363 events in the sideband region. 
The latter, together with the effective ratio of areas of the signal region to the sideband region, 
implies an expected number of background events in the signal region of 6.24 it 0.33 it 0.36. The 
quoted systematic uncertainty comes from the background shape variation discussed previously. 
Figure 4 shows a projection of the data on to the mss axis after requiring \AE\ < 25 MeV. 

We use a detailed Monte Carlo simulation of the BABAR detector to determine the efficiency for 



reconstructing the signal This, together with the total number of BB pairs produced during 
data collection, allows us to determine a preliminary branching fraction for B^ D*^ D*~to be 
5(5° ^ D*+D*-) = (8.0 ± 1.6 ± 1.2) X 10"! 
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Figure 3: Distribution of events in the AE vs. mES plane. The small box indicates the signal 
region, while the sideband region is everything outside the larger box. 



(N 

o 



in 

^— < 

> 



10 







BaBar 



5.2 



5.22 



5.24 



5.26 



5.28 



(GeV/cO 



Figure 4: Distribution of events in tties plane with a cut of jAE'l < 25 MeV applied. The curve 
represents a fit to the distribution of the sum of a Gaussian to model the signal and an ARGUS 
function to model the background shape. 



The dominant systematic uncertainty in this measurement comes from our level of understand- 
ing of the charged particle tracking efficiency (9.4%). The high charged particle multiplicity in 
this decay mode makes this measurement particularly sensitive to tracking efficiency. Uncertainties 
were assigned on a per track basis for vr, K and slow tt, and were added linearly due to large 
correlations. The imprecisely known partial- wave content of the S° D*+D*-&nal state IS an- 
other source of systematic uncertainty (6.6%). This was estimated by calculating the change in 
the reconstruction efficiency for different final angular states in Monte Carlo. Other significant 
systematic uncertainties arise due to the uncertainties on the D*^ , and branching fractions 
(5.6%) and the differences in mass resolutions between Monte Carlo and data (4.1%). The total 
systematic uncertainty from all sources is 14.5%. 



3 Summary 

Using about 23M BB events, we have observed several hundred completely reconstructed B — >■ 
D^*^D^*^ K decays. The following preliminary branching fractions have been measured: B{B^ 
D*-D^K+) = (2.8 ± 0.7 ± 0.5) x lO'^ and B{B^ D*-D*^K+) = (6.8 ± 1.7 ± 1.7) x 10-^ in 
good agreement with the CLEO measurements B{B^ D*-D°K+) = (4.5l?;9 ± 0.8) x 10"^ and 
B{B'^ D*-D*^K+) = (13.0+^ j ± 3.6) x 10^^ |l|]. We have observed an excess of 8.2 ± 3.5 



events over a background of 1.7 events in the color-suppressed decay mode B^ — > D*^D* K~^, 



where no significant number of candidates has been previously seen. The corresponding prehminary 
branching fraction is measured to be 13{B~^ D*~ D*^ K^) = (3.4 ± 1.6 ± 0.9) x 10"^. Finahy, 
several candidates have also been observed in the CP eigenstate D^*^^ D^*^^ Kg. This study 

confirms that the transitions b — > ccs can proceed through the decays B D^*^D^*^ K. We have 
observed a signal of 31.8 it 6.2 it 0.4 events in the decay B^ — > D*^D*^ . We measure a preliminary 
branching ratio to be B{B^ D*^D*~) = (8.0 ± 1.6 ± 1.2) x 10~^ in good agreement with the 
CLEO measurements B{B^ D*+D*-) = (9.9t|^ ± 1.2) x lO""^ |T|]. 
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